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Abstract 

In this article effect of radiation preheating of unburnt mixture by prop- 
agating deflagration front is studied from the viewpoint of its ability to form 
a promoting temperature gradient and trigger transition to detonation. First, 
we investigate the effect of radiation preheating of the unburnt mixtures, 
when they are traveling through the wrinkles on the flame surface, in or- 
der to estimate a possibility of significant temperature rise. Subsequently, 
numerical simulations of a simplified mathematical model are carried out. 
They demonstrate plausibility of the proposed mechanism of the deflagra- 
tion to detonation transition. 

Keywords: Premixed combustion, transition to detonation, DDT, radia- 
tion heat transfer. 



1 Introduction 

If a large amount of accidentally released gaseous fuel is within the limits of 
fiammability and ignition takes place, a chance for transition of the deflagration 
front into the detonation wave (DDT) appears. Actualisation of the DDT depends 
on a variety of circumstances and may develop in a number of physical scenarios. 
Among all of these scenarios, the mechanisms not involving interference from 
effects, substances and objects external to the released fuel, i.e. incident shock 
waves, jets, particular geometry of nearby walls and objects, pose special interest. 
This is virtually the only possible mechanism of DDT in open atmosphere, and 
thus is especially important for safety in large storage and processing petrochem- 
ical plants and pipelines, see e.g. 0]. 
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Research in explosion safety usually assumes that principles of safety by de- 
sign rule out powerful sources of ignition in explosion prone industrial environ- 
ment leaving chances to more delicate mechanisms of the DDT. For example ob- 
stacles placed in the way of the deflagration front or incident shock waves, see 
e.g. ElEl. Alternatively, they can be self-initiated. 

Probably the first successful attempt to explain possible mechanism of self- 
initiated DDT was the idea of the promoting temperature gradient. According 
to JH, temperature gradient in the reacting media results in a gradient of induc- 
tion times of chemical reactions. If the gradient of the latter matches the speed of 
propagation of the reaction front, then the synergy of the two processes may trans- 
form the flame into a detonation wave. This scenario was carefully investigated 
in a number of later studies, see e.g. [5 ], where an idealized medium with simple, 
rate-sensitive kinetics and an initially prescribed linear gradient of temperature 
was studied exhaustively. 

Once it was understood that temperature nonuniformities may induce transi- 
tion to detonation, an investigation of possible physical mechanisms leading to 
such uniformities was undertaken. Amongst theoretical achievements in the area 
we would mention the concept of the equivalent drag force, which works to cre- 
ate temperature nonuniformities ahead of the flame front J6]|. Numerical simu- 
lations demonstrating importanve of the effect of viscous friction in the bound- 
ary layer were carried out in relatively narrow channels. Although, further 
work [0 indicates that the effect remains relevant in quite realistic cases too. 

Another approach is based on the collective thermal effect of shock waves 
generated in front of a very fast flame by the flame itself JU. The theory looks 
feasible for flames, which are naturally very fast and, in addition, experience im- 
mense exponential acceleration. Both computational and experimental examples 
given in ||9l demonstrate possibility of such exponential acceleration of flames in 
tubes up to supersonic propagation speeds with Mach numbers well in excess of 
one. At such speeds supersonic deflagration fronts are able to generate down- 
stream shocks strong enough to change thermal and reactive properties of the fuel 
cardinally. It is not clear how to extend these ideas to unconfined outward propa- 
gating flames. 

1.1 Acceleration of expanding flames 

It was discovered in ifTOll . that expanding flames do accelerate. This accelera- 
tion was linked to appearance of wrinkles, as wrinkles increase flame surface 
area allowing to burn larger quantities of fuel and thus expand faster. Theoretical 
and numerical analysis, see e.g. IfTTT - TOll . revealed physical mechanism and basic 
characteristics of the acceleration. It was established that flame fronts become ex- 
tremely sensitive to the upstream velocity perturbations as their size grows. From 
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the physical point of view, the sensitivity of cellular flames to perturbations of 
the upstream flow velocity is explained by the possibility of a pseudoresonant 
interaction between them lfl4ll . The pseudoresonant interaction amplifies certain 
upstream velocity perturbations and transforms them into pairs of coupled vor- 
texes of size proportional to the the thermal flame thickness, which are located 
just behind the flame front. These pairs of vortexes suck the unburnt gas into the 
areas occupied by the products of combustion and form wrinkles which signifi- 
cantly increase the bulk burning rate and hence the averaged flame propagation 
speed. 

The power law of growth of the averaged flame radius f(t) oc (t — i*)^ for large 
enough times t > t*, in accordance with the fractal analysis IPT51I . was confirmed. 
Values of the power depend on the dimension of the problem and equal 3/2 for 
three dimensional flames and 5/4 for two-dimensional ones. In expanding flames, 
the pseudoresonant interaction between the flame and the upstream flow velocity 
perturbations is weakened by stretch. As a result set up of the acceleration at i* 
can be delayed significantly. On the other hand, permanent growth of the size 
of the flame results in its permanent acceleration unlike the non-stretched planar 
flames, which have a finite limit on their propagation speed [[TBI . Using a realistic 
set of dimensional parameters, e.g. planar flame speed relative to the burnt gases 
Ub = 0.5 m/s and thermal diffusivity D t h = 2.5 x 1CT 5 m 2 /s, our findings for t* in 
dimensional terms can be put as f(t*) ~ 20 m. 

Theoretically, continuous acceleration of the flame should lead to supersonic 
speeds and further to DDT. However, the rate of expansion of spherical flames 
f(t) oc (t — t*) 3//2 , is so slow that even an optimistic estimation, e.g. assuming 
initial flame expansion rate of about v « 66 m/s for flame of averaged radius r ~ 
20 m, gives flame radius f of order of ro(v/v ) 3 « 2.5 km for which supersonic 
speeds v ~ 330 m/s can be reached. This means that if self-initiated DDT of 
expanding flames takes place at all, it should be assisted, and even more likely 
dominated, by an effect other than flame acceleration. 

1.2 Radiation preheating 

Upon realizing that flame acceleration due to flame cellularization alone cannot 
lead to self-initiated DDT, appearance of a promoting temperature gradient in 
front of expanding flames looks the only factor which can. The question is how 
such temperature gradient can be generated in front of the expanding flame? Con- 
sideration of time scales specific to heat transfer via conduction and radiation, 
shows that the most likely heat transfer mode which is able to modify the tem- 
perature gradient is radiation. Normally, gaseous fuels do not absorb radiative 
heat significantly, however even minute quantities of water vapour can change 
this behaviour dramatically. Moreover, fine aerosol and dust fuels often are good 
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radiation heat absorbents. Initiation of detonation by ultraviolet radiation was 
demonstrated in IfToll . where the radiation was an external factor and its direct 
action was photochemical rather than the thermal one. 

Calculations show that the depth of the wrinkles on two-dimensional circular 
flame surface reaches up to 1/10-th of the averaged flame radius. Comparison 
with the three-dimensional spherical flames shows that the latter expand and grow 
cellular structures even faster, see Fig. Q]from lfT2ll . Their cellular structures are 
better developed and wrinkles are even deeper. This suggests that some unburnt 
fluid particles may spend significant time inside the wrinkle before reaching the 
flame front itself. Thus, preheating of the unburnt mixture, traveling inside the 
wrinkles, with the heat radiation from surrounding flame surfaces might be im- 
portant. 

xio 5 
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Figure 1: Development of a fractal structure on surface of circular flame (left). 
Comparison of the circular flame and the equatorial cross-section of the three- 
dimensional spherical flame at nearly the same time instances (right). 

In this work we carry out tentative estimations of a possibility of the effect of 
radiation preheating to trigger the DDT. In Section [2] we consider a simple ana- 
lytical model of the phenomena confirming a possibility of significant growth of 
temperature due to radiation in deep enough wrinkles. Further, in Section [3j a 
Navier-Stokes based numerical model is presented. Results of a set of computer 
simulations of the phenomena in question are discussed in Section 0] indicating 
that the DDT in combustible mixtures with essential heat emissivity and absorp- 
tion is possible. Physical mechanism of such heat radiation triggered self-initiated 
DDT's is explained too. 
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2 Preliminary estimations 



2.1 Heat balance equation 

Let us assume that a combustible fluid particle P is moving inside a flame wrin- 
kle and estimate its temperature growth because of heat radiated from hot flame. 
The estimation will be by order of amplitude, similar to that one outlined in ifTTl . 
For the sake of simplicity the products of combustion are assumed at the same 
uniform temperature T w and particle P is not very close to the flame yet. Under 
this conditions, the only heat which matters is the heat from the flame absorbed 
by the particle. All other elements of the radiating system are in nearly thermody- 
namical equilibrium and radiation heat exchange between them can be neglected. 
Furthermore, the wrinkle will be modelled as a conical cavern of depth L and of 
the apex semi-angle a, see Fig. 13 Under assumptions in question, variations of 
the geometry of the wrinkle result just in a coefficient of order one and are not 
essential. 




Figure 2: Model of a wrinkle on the surface of an expanding flame. 



Every element W of area dA w of the cone surface radiates heat according 
to the Stefan-Boltzmann law with the coefficient of emissivity ew- So that heat 

dE 

= Ew&TwdA w . The heat is transferred from W 



radiated per unit time is — — 

F dt 

to P, which is at a distance /, dissipating according to the exponential law with 



w 
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the dissipation factor \- It means that heat radiated from point W per unit time 
will be weakened by e~ PaXl times upon travelling through distance I. Here p a 
is the density of the absorbing species, which might be quite different from the 
averaged density p of the combustible mixture. 

The heat flux incident on the fluid particle P of cross-section area dAp is the 
areal fraction dA P / (2nl 2 ) of the radiation energy which reaches P after being 

~"" xl sw^T^dAw- Hence, the fraction of 



dE 

radiated from W is — — 

dt 



2ixl 2 



this radiation heat which will dissipate within the small fluid particle P of thick- 

P (p a xdl) e~ paXl (ew&Tw) dA w . Denoting volume of 



... dE 
ness dl is — — 

dt 



2tx\ 2 



the fluid particle P as dV = dApdl, the total rate of heat gain of the fluid particle 
Pis 



dE P _ ZwVpaXT w f j-2-paxUA 
dVdt ~ 27T J 1 6 dAW - 



A 



w 



Noticing that the x, ^-averaged I 2 = (£— x) 2 + (£ — L) 2 (tana) 2 is proportional 



L — £ 

to L 2 , dA w = d£d(fi, and that the integrand of the above integral is positive, 

cos a 

the integral itself can be approximated as follows 



dE P ^ e w ap aX T^e-^ (]> 
dVdt 2c g cos a 

where c g is a "geometry correction coefficient" of order 0(1). 

As temperature of the fluid particle grows due to radiation heating, inten- 
sity of chemical reaction starts to grow too. We estimate the rate of this heat 
release assuming an overall exothermic uni-molecular irreversible chemical reac- 
tion A — >■ B. For example, A is a mixture of CH A and air and B is a mixture 
of air, C0 2 , and gaseous H 2 0. If the reaction produces energy Q per mole of 
reacted species A, then the heat release because of chemical reactions in the fluid 

particle per unit volume per unit time is p ^ = Qp p k a e^^ , where cp is the 

specific heat capacity at constant pressure, v is the order of the reaction, k a is the 
pre-exponential factor, and E a is the activation energy. Assuming that the fluid 
particle is moving toward the flame front with velocity u n , we can replace time t 
with x = u n t: 

dT ^ ... Eg 

cppu n — = Qp k a e kt ; (2) 
dx 

Combining heats produced by radiative heat flux from the flame surface © 
and chemical reactions ©: 

dQ TaO 

— « g c ei+e +q r , (3) 
dC 
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Qpr- 1 k a Le- T ° e w aT^ PaX Le-P^ L V^ E 

where q c = , q r = , T a = ——, 6 = 

cpu n ±Q c g pcpu n l o cos a HIq 

rp ° . and C = x/L. 
2.2 Approximate solution 

In the early stages of heating, fluid particle temperature does not differ from T 

( T a Q 

significantly and the Arrhenius exponent can be approximated as exp f - — — 
e T a G j^en^ correspondingly approximated heat balance equation ([3]) is inte- 

1 + ^l) e -*rT a C _ ^ 

q r J q r _ 



grated exactly @(£) ~ ——hi 

J- n 



yielding the solution 



blowup distance for £ = x/L equal 



Cw— r ln(l + ^). (4) 

q r T a V QcJ 

It is interesting to compare the latter to the blowup distance Q ~ — f° r Qr = 0. 

Considering parameters typical to combustion of hydrocarbons in air T a = 80, 
Q/(cpT ) = 8, u n = 0.5 m/s, fc a = 5 x 10 12 s^ 1 , v — 1, one obtains g c w 
2 x 10 _22 L. Taking into account that function f(x) = xe~ x reaches its maximum 
equal e _1 for x = 1, we conclude that q r is maximal for p a xL = c g and this 
maximum is 

qr, max ^3/2 X m " v ^ 

ec 9 cos a pcpUn 1 ^ 

Assuming that the first factor in © is about 1/2, p = 1 kg/m 3 and T w = 1500 K, 
one obtains from © that w 0.9 for L > lm. Numerical integration of © 
shows that significant growth of temperatures can be reached for even smaller 
flames. According to our earlier analysis, wrinkles of depth 1 m can appear on 
flames of radius r c ^ 10 m. 

~ —1/2 

Our choice of p a xL — °g ma y artificial. However, it is not, because 
the value of the product is effectively controlled by the density p a of the heat 
absorbing admixture and by the wrinkle depth L. The former, is at the disposal of 
the experimentalist. For some particular species values of p a suitable for efficient 
combustion might be limited and the limit may lie below the requirements of the 
efficient radiation heat absorption for given \ an d L. Even so, required values of 
p a can be always lowered below this limit for deeper wrinkles, i.e. for larger L. 

These simplified estimations were carried out for some "typical" values of 
combustion and heat transfer parameters and for an optimal radiation heat ab- 
sorption. Thus, they just show that radiation preheating of fuel in front of the 
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deflagration flame might be significant in principle. Such temperature rise alone 
not necessarily leads to the DDT yet and more detailed numerical simulations are 
still required in order to validate the hypothesis. Furthermore, mixtures with such 
a combination of the parameters in question might be difficult to find in practice 
even if theoretical numerical simulations confirm that radiation preheating can 
cause the DDT. 

As we have mentioned in Section 11.21 two heat transfer mechanisms poten- 
tially able to produce promoting temperature gradient are heat conduction and 
radiation. According to ©, the time scale of the radiation preheating over dis- 
tances of order l a = (paX) 1 1S °f order r r = ^ Cp ^ 3 . In contrast, the time scale 

of heat transfer by conduction over the same distance is r c = l 2 a / n, where k is 
the thermal diffusivity. By representing k oc 8 t hU n in terms of the thermal flame 
thickness 8 t h and laminar burning velocity u n , the ratio of the time scales takes 
the form 

T r _ pC P U n 8 t h _i 

r ~ aT 3 X / £w ' 

1 c u - L W a 

For gases value of pep is typically about 10 3 J/(m 3 xK) and u n is of order 
lm/s. Then, value of the Stefan-Boltzmann constant a w 5.67 x 10~ 8 kg/(s 3 xK 4 ) 
and reasonable flame temperatures of at least T w = 1800 K ensure that the first 
factor in © does not exceed 3. Hence, for typical gaseous fuels 

Tr 3 8 t h f£S 

— = — x — . (6) 

Thickness of the preheated layer l T , which renders optimal temperature gradi- 
ent, exceeds 8 t h by a few orders of magnitude Q. Thus, for reasonably radiating 
flames, i.e. ew > 0.03, and if we manage to adjust p a to make the dissipation 
length l a match It, radiation heat transfer will dominate conductivity. 

Relationship between the leading time scales is getting more complicated for 
multi-component and especially for multiphase fuels. For example, if radiation 
absorbing and reacting species are different, then the speed of heat transfer be- 
tween them should be taken into account too. Further complications are linked 
to the time scales of evaporation and pyrolysis processes, which are relevant to 
burning of mist and dust mixtures. 



3 Mathematical model and its numerical implemen- 
tation 

Numerical simulations of dynamics of deflagration fronts in presence of inten- 
sive radiation preheating was studied in the case of planar flames propagating in 
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tubes open from one end rather than for expanding spherical flames in the open. 
This geometrical simplification made it possible to carry out investigations with 
relatively limited resources. 

The approach used Navier-Stokes description of compressible fluid and em- 
ployed drastically simplified chemistry and radiation models. The governing equa- 
tions can be written as follows 

ov + ar(vo + <W0 = 0.7*0 (vq + ag^jv) + 

dt dx\ 0x2 dxi dxi 



where V 



p, pu x , pu 2l ^ + ?^fp {u\ + u 2 2 ) , pY 



Formulas for convective 



J~(V), Qiy), and dissipative F^iV), Qw(V) fluxes can be found elsewhere. 
The source term Q(V) is combined of contributions from chemistry and radiation. 
Tube diameter d and laminar burning speed relevant to the burnt gases u b were 

Ufa 

used as scales, so that the Mach number M — — , where a is the speed of sound 

a 

in the fuel at initial temperature T . Other nondimensional governing parameters 

used were the Reynolds number Be = P° Ub ^ ^ T a = and Q = . Prandtl 

P RTq ' c v T 

Pr = and Lewis Le = numbers were set to one. 

K CpDp 

Chemistry was modelled with a single irreversible A B reaction of Ar- 
rhenius type similar to the assumptions in Section [27TI The radiation model was 
reduced to the heat flux driven by temperature according to the Stefan-Boltzmann 
law. In addition, emission by cold gas and gas too far away from the flame as well 
as absorption by hot gas and gas too far away from the flame were not taken into 
account, because gases in those areas are in nearly thermodynamical equilibrium 
from the view point of radiation. In spite of its roughness, the model governs 
interaction of all relevant energy fluxes reasonably well. 

The governing equations of the model were solved with a high order upwind 
shock capturing finite difference scheme without any fractional steps. Numeri- 
cal approximation is implicit in the stream-wise coordinate, and is explicit in the 
other one, similar to ffT8l . providing satisfactory compromise between limitations 
of numerical stability and efficient parallelization. Computational domain was 
moving together with the flame front and absorbing boundary conditions were 
implemented on its up- and downstream boundaries. The resulting code is reason- 
ably efficient for Mach numbers greater than 0.01, Reynolds numbers less than 
1000, and for the nondimensional activation energy T a less than 50. 
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4 Computer simulations and discussion 



4.1 Numerical experiments 

In our numerical simulations self-initiation of DDT never took place in smooth 
open pipes without obstacles and absorption of radiation heat. On the other hand, 
significant absorption of radiation heat was able to alter behavior of the flame 
front dramatically. It was found that self-initiation of DDT in radiation heat ab- 
sorbing media progresses via formation of an optimal temperature gradient on the 
sides of wrinkles on the flame surface. Successful DDT was resulting from inter- 
action ("collision") of two flame segments accelerating through such promoting 
temperature gradients as illustrated in Fig [3l Zone of intensive chemical reaction 
is indicated with a solid line in there, while the boundary of the layer of unburnt 
gas at essentially elevated temperature is shown with a dashed line. Normals to the 
flame surface indicate presence of parts of the preheated layer in which structure 
of the temperature gradient is the most advantageous for the deflagration front to 
transform into a detonation wave. Appearance of such a thick preheated layer of 
the unburnt gas can be seen in the second row in Fig. [5] as well. 




Figure 3: Formation of an optimal temperature gradient by radiation preheating. 

Sufficient amount of fuel available between two opposing flame segments of 
the wrinkle is critical for success of the DDT. In addition, the effect is sensitive to 
the shape of the flame front, which in turn depends on a set of governing param- 
eters of the problem. For example, chemical reactions of orders higher than one 
widen range of suitable flames able to experience the DDT significantly, confirm- 
ing earlier observations of other researchers, see e.g. |[T9ll . Note, heat radiation 
from the flame not only heats up the unburnt fuel in front, but cools down the 
flame itself too, which levels out the gradient even faster. 

In contrast to expectations based on common sense, success of the self-initiated 
DDT depends on the depth of the flame wrinkle in a lesser degree, which makes 
our estimation of critical flame radius of order r c « 10m in Section [2721 less rele- 
vant. The self-initiated DDT based on the mechanism just described might work 
for much smaller flames. An asymptotic estimation of the critical radius of out- 
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ward propagating unconfined flames for this new self-initiated DDT mechanism 
is possible, but would deserve a separate investigation. 

As absorption rate grows, flame dynamics evolves through a sequence of dis- 
tinctive regimes. First, the flame is not affected by the radiation at all. Rise of 
the absorption rate leads to formation of promoting temperature gradients able to 
generate powerful compression waves. However, these temperature gradients are 
insufficient to allow for proper acceleration of deflagration fronts yet and gener- 
ated compression waves are not strong enough to maintain supersonic combustion 
behind them. As a result, the steady flame propagation is altered by a sequence of 
failed detonations as illustrated in Fig. 0] 




5.5 6 6.5 7 




5.5 6 6.5 7 



Figure 4: Example of evolution of pressure (left) and reaction rate (right) fields in 
the failed self-initiation of DDT. 



Eventually, when absorption rate is high enough, stable self-initiated DDT 
takes place GUI as this is illustrated in Fig [5J The first row depicts the flame just 
before the DDT. The next row illustrates formation of the preheated layer of the 
unburnt gas in front of the flame. Soon after events begin to unfold really fast. 
The third row shows the moment when a "micro-explosion" inside the wrinkle 
just took place. Generated powerful compression wave kick starts the detonation 
front shown in the following rows. Reminiscences of the shock waves brought 
forth by the opposing flame segments continue to interfere with the downstream 
propagating detonation front for some time as can be seen in the forth row of Fig. 
\5\ However, they decay eventually and the detonation wave settles gradually to a 
perfect ZND detonation structure exemplified in the fifth row. 

Effect of radiation heat transfer at very high radiation absorption rate was not 
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Figure 5: Example of evolution of pressure (left) and reaction rate (right) fields in 
the self-initiation of DDT. 

studied yet, but is believed to be similar to the case of high thermal conductiv- 
ity. Hence, dynamics of the self-initiated DDT at very high radiation absorption 
rate is expected to be similar to the case when optimal temperature gradient is 
formed by means of heat conductivity rather than by the radiation heat transfer 
©. It is interesting to note that the detonation initiations presented in [fT9ll were 
obtained for rather high thermal conductivity and that the range of flame wrinkles, 
able to generate a detonation, widens as heat conductivity grows. All this can be 
explained in terms of formation of an optimal temperature gradient by means of 
conduction. It is also known that at high thermal conductivity flames lose abil- 
ity to form properly shaped wrinkles and an artificial corrugated wall was needed 
in lfl9l to create them. Introduction of strong fuel diffusion, i.e. consideration of 
small Lewis numbers, see e.g. ETI . might fix the problem too. 

4.2 Critical analysis 

Practical realization of the self-initiated DDT mechanism in question is expected 
for mixtures containing water vapour, fine liquid sprays and dusts because vari- 
ation of concentration of the admixtures affects values of the dissipation factor 
PaX m a very wide range. Many experimentalists report that moist atmospheres 
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increase chances of DDT significantly, see e.g. E2ll . However, in practice, choice 
of mixtures with required properties might be limited. Moreover, mixtures with 
specific characteristics, which can be easily considered theoretically, may simply 
not exist in real at all. Hence, validation of practical usefulness of the proposed 
mechanism of DDT requires development of realistic chemical and heat radia- 
tion/absorption models. These models should be coupled, i.e. changes in the con- 
centration of water vapours should affect not only rates of heat radiation and ab- 
sorption in the mixture, but rates of chemical reactions as well. High enough con- 
centration of water vapours will result in almost as high heat radiation/absorption 
rate as required, but at certain stage their excess will begin to slow down or even 
completely extinguish the flame. 

On the other hand, once quite detailed chemistry is going to be considered, the 
heat radiation/absorption model has to be refined too. As concentration of water 
vapour changes, new intermediate species with special radiation properties may 
appear. Therefore, the heat radiation/absorption model should take into account 
the spectral structure of the radiation fluxes rather than their total energy alone. 
Furthermore, if heat absorbing and reacting species are not the same, then heat 
exchange between them should be taken into account too. 

Relatively low Reynolds numbers and activation energies were used in our 
numerical simulations. Character of changes as the Reynolds number grows does 
not suggest any problems up to transition to turbulence. The latter looks just a 
technical challenge rather than the undoer of the proposed DDT scenario. Changes 
in the character of flame dynamics when activation energy grows suggest just 
further shrinking of the range of absorption rate for which detonation repeatedly 
fails. The tendency explains lack of experimental observations of the regime of 
successively failing detonations. 

Eventually, this investigation is based on a simplified geometrical model of the 
phenomena. The original problem of self-initiated DDT is a three-dimensional 
one and is in the open atmosphere, whereas only a planar two-dimensional com- 
bustion in a pipe was considered in this paper. 

5 Conclusions 

Possible temperature rise of gas fuel in front of propagating flame as a result of ra- 
diation heat transfer was estimated. The estimation demonstrated that temperature 
rise of fuel approaching flame inside a wrinkle can match the adiabatic burning 
temperature for relatively small flames of radii just 10 m as soon as the gas mix- 
tures involved are able to emit and absorb radiation heat at high enough rate. 

A variety of numerical experiments on flame dynamics in pipes open from one 
end was carried out as well. As absorption rate grows, flame dynamics evolves 
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through a sequence of regimes. Initially flame is not affected by the radiation. At 
certain rate of radiation heat absorption flame propagation is altered by a sequence 
of failed detonations. Eventually, stable self-initiated DDT takes place. 

Undertaken simulations demonstrate plausibility of the idea of the radiation 
preheating as the principal effect in the self-initiated DDT . However, the self- 
initiated DDT is a very delicate phenomena. More accurate and detailed numerical 
simulations are needed in order to clarify issues highlighted in this work. 
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